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Abstract 

We observe a hyperfine anomaly in the measurement of the hyperfine splitting 
of the 6S1/2 excited level in rubidium. We perform two step spectroscopy using 
the 5Si/ 2 — > 5P]/ 2 —* 6S1/2 excitation sequence. We measure the splitting of 
the 6S1/2 level and obtain for the magnetic dipole constants of 85 Rb and 87 Rb 
A = 239.18(4) MHz and A = 807.66(8) MHz, respectively. The hyperfine anomaly 
difference of 87^85 = —0.0036(2) comes from the Bohr Weisskopf effect: a correction 
to the point interaction between the finite nuclear magnetization and the electrons, 
and agrees with that obtained in the 5jS]/2 ground state. 
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Measurements of hyperfine splittings provide a low energy approach to study 
the magnetic moment distribution in the nucleus. The hyperfine splitting is 
due mainly to the coupling between the magnetization of the nucleus with 
the magnetic field created by the electrons. To a very good approximation 
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the coupling is taken to be point-like. However, high precision measurements 
of the splitting can reveal deviations from the point interaction. These de- 
viations, also called hyperfine anomalies, come from considering how finite 
magnetic and charge distributions affect the interaction. Bohr and Weisskopf 
[lj first discussed the finite magnetization effect in the anomaly. The modified 
charge potential that the valence electron sees as it gets closer to the nucleus, 
the Breit-Crawford-Rosenthal-Schawlow [2113] effect, is the other source for a 
hyperfine anomaly. 

We show in this letter a measurement of the excited state hyperfine split- 
tings in two isotopes precise enough to extract a difference in the hyperfine 
anomalies. The hyperfine anomaly differences of the ground state in francium 
were measured by our group [I] by comparing the ratio of the hyperfine split- 
tings of the 7P\/2 (which has a small anomaly) and the 7S\/2 (which has the 
large anomaly) [5]. These measurements represent excellent benchmarks to 
test state-of-the-art ab initio calculations of the electronic and nuclear wave 
functions for future parity non-conservation measurements [BUT] . They give the 
possibility to look at how the systematic addition of a neutron modifies the 
electronic wavefunction in the nucleus, and how the nucleus rearranges itself. 
As a new generation of proposed parity violation experiments in atoms starts, 
PI51I9] . it is important to understand the nuclear structure limiting factors for 
the extraction of weak coupling parameters from the measurements [TU] . This 
letter presents a new experimental approach to extract and study such effects 
in chains of isotopes. 

The hyperfine shift for these levels is given by Ehf = A(F(F + 1) — 1(1 + 
1) — J (J + l))/2 where A is the magnetic dipole constant, F is the total 
angular momentum, I is the nuclear spin, and J is the total electronic angular 
momentum. Derivations of A assume that the nucleus is a point particle with 
magnetic moment fijy and charge Z. However, this is not the case. The nucleus 
has a finite extension and structure. The value of the magnetic dipole constant 
has to be modified to include the effect of the charge and magnetic distribution 
on the electronic wave function. We can write, following Ref. [TTJ A for an 
extended (ext) nucleus as a function of the point value (point) (Eq. [[]): 



Apoint = —J^gi^N^B\^P(0)\ 2 fR, 
A ex t = A point (l + €bCRS){1 + tBw), (1) 

where ip(0) is the electronic wave function evaluated at the center of the nu- 
cleus, fiB is the Bohr magneton, fi N is the nuclear magneton, gj is the nuclear 
g-factor, and represents the relativistic enhancement. The last two terms 
of this expression modify the hyperfine interaction to account for an extended 
nucleus. The first of them (escRs), the Breit-Crawford-Rosenthal-Schawlow 
(BCRS) correction, accounts for the fact that the valence electron sees a mod- 
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ified potential as it gets closer to the nucleus. The second one (esw), the Bohr 
Weiskopf (BW) effect, focuses on the magnetic distribution of the nucleus. 



The fractional difference in the mean charge radius between the two isotopes of 
interest ( 85 Rb and 87 Rb) is less than 10~ 3 and justifies neglecting the BCRS 
correction [T2"1[T3"] . 87 Rb has a closed-neutron shell that makes the nuclear 
charge distribution insensitive to the addition or subtraction of neutrons |12j . 
A calculation of €bcrs P3] gives a value of 1% compared to the point nucleus; 
however, since we are interested in comparing the two anomalies we find, 
using the values of the charge radius of Ref. [13], the difference between the 
corrections of both isotopes is less than 10~ 4 . 

The BW correction describes the modification of the hyperfine interaction due 
to a finite distribution of nuclear magnetization. The addition of more neutrons 
to the nucleus yields a different spatial magnetization. The magnitude of this 
correction for an individual nucleus is smaller than the BCRS correction but 
depends heavily on the spin and orbital angular momentum of the nucleus and 
not just on its radius. For the particular case of 85 Rb and 87 Rb the addition of 
two neutrons changes the value of the valence proton from an fen orbital to a 
P3/2 orbital. The lighter of the isotopes, 85 Rb, has the spin angular momentum 
of the valence proton anti-aligned with its orbital angular momentum while 
87 Rb adds them both together. The two neutron holes deform the nucleus 
very slightly and change the order of the orbitals that are almost degenerate. 
The difference of the hyperfine anomaly contributions 87^85 = e|V — e §V * s 
dominated by the B W contribution as we take the BCRS correction as equal 
for both isotopes. 

Since the magnetic moments of the nuclei are well known [15] we can use a 
ratio to extract the Bohr-Weisskopf effect difference: 

A 87 a 85 

— ~l+ s7 5 85 . (2) 
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Because the BW corrections are small, we require a method to measure the 
hyperfine splitting with good precision. We perform two-photon absorption 
spectroscopy in a 30 cm long 2.5 cm diameter Rb cell with natural isotopic 
abundances to measure the hyperfine splitting of the 6Si/ 2 level using the 5Pi/ 2 
as an intermediate level (see Fig. [1] for the block diagram of the experimental 
setup and the atomic energy levels involved). The cell has no buffer gas and 
is in a controlled magnetic environment. 

We monitor the change in absorption of a laser resonant with the transition 
to the 5Pi/2 level at 795 nm as we scan the 1.3 /im laser over the hyperfine 
structure of the 65*1/2 level. We modulate the 1.3 /im laser at a frequency 
close to half of the hyperfine separation to produce sidebands well separated 
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A- (Modulator} — 1 .3 urn laser 




Fig. 1. Schematic of the experimental apparatus. The inset shows the energy levels 
relevant to our experiment (not drawn to scale) with quantum numbers correspond- 
ing to 85 Rb ( 87 Rb). 



in energy. This way we only need to scan the laser a small range around 
the midpoint between the two hyperfine levels of the 65*1/2 level such that 
the sidebands scan over both hyperfine levels. The sideband signals appear 
in the absorption profile and work as an in situ scale. We record absorption 
profiles at different values for the modulation to find, with the help of linear 
regression plots, the point where the sidebands cross the midpoint of the two 
main peaks. Typical fits give correlation coefficients that differ from unity by 
at most 2 x 1CT 4 . 



A Titanium Sapphire laser with a linewidth better than 500 kHz tuned to 
the Rb Di line at 795 nm excites the first step of the transition. A Pound- 
Drever-Hall setup on an independent Rb cell locks the 795 nm laser to the 
F=l(2)— >F'=2(3) transition in 87 Rb ( 85 Rb). A grating narrowed diode laser at 
1.324 fim with a linewidth better than 500 kHz excites the second transition. 
We scan the frequency of the 1.3 /mi laser using a triangular shaped voltage 
ramp from a synthesized function generator at 4 Hz applied to the piezo control 
of the grating and monitor its frequency with a wavemeter. A fiber-coupled 
semiconductor optical amplifier increases the power of the 1.3 fim laser before 
it goes to a large bandwidth (« 10 GHz) Electro-Optic Modulator (EOM). A 
signal generator modulates this EOM in the range between 100 and 900 MHz 
to imprint the sidebands necessary for the measurement. 
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Before going through the rubidium vapor cell the laser beams are circularly 
polarized by A/4 waveplates. A counter propagating 1.3 /iin laser beam with 
a power of up to 4 mW overlaps one of the 795 nm beams. The relative power 
on the sidebands is less than a tenth of the total and its waist is 1 mm. The 
lasers overlap to a precision of better than 1 mm along 75 cm. The cell resides 
in the center of a 500-turn solenoid (7.4 Gauss/A) contained inside a three 
layered magnetic shield. The solenoid is 70 cm long and has a diameter of 11.5 
cm. 

A thick glass plate splits the 795 nm laser beam into two co-propagating 
beams. The power of each beam is approximately 10 /xW with a waist of 1 
mm. We find that doing the two step excitation in either a (cr + , a~) or (a~, <r + ) 
polarization sequence for the 795 nm and 1.3 fim lasers respectively, increases 
the probability of the transition to the 65*1/2 level, and allows us to check for 
optical pumping and Zeeman shift effects. We place the Rb cell in a uniform 
magnetic field convenient for the Amp = ±1 transitions. The addition of the 
magnetic field has the two-fold advantage of providing a quantization axis 
and a tool to probe systematic effects. We measure the hyperfine splitting for 
different values, including reversals, of the magnetic field and polarization. We 
can then extract by interpolation the value of the hyperfine splitting at zero 
magnetic field. 




Fig. 2. Experimental traces that illustrate sideband crossing for Rb. The intersec- 
tion of the two projected lines marks half the hyperfine separation. 

After the glass cell an independent photodiode detects each 795 nm laser 
beam. The output of both detectors goes to a differential amplifier to reduce 
common noise. A digital storage oscilloscope records the signal and averages 
for about three minutes. The stability of the 1.3 /iin laser center frequency 
over this time is better than 0.3 MHz as measured by a Fabry Perot cavity. 
The optical attenuation at line center (Di line) is 0.4 for 85 Rb. Fig. [2] presents 
an example of the oscilloscope traces for different sideband frequencies. The 
projection shows the positions of the peaks joined by the linear interpolation, 
the error in the position is much smaller than the dots. We obtain signal to 
noise ratios in excess of 100 for the sideband resonances. 
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We analyze how the following parameters influence the hyperfine separation: 
the peak shape model for the non-linear fit to obtain the separation of the cen- 
ters of the lines, scan width and repetition rate of the 1.3 /xm laser, power of 
the 795 nm and 1.3 /im lasers. We look for optical pumping effects, magnetic 
field effects, and temperature, though the vapor density, that can influence 
the hyperfine separation. We have only been able to put upper bounds on 
the systematic shifts within our resolution and we determine that statisti- 
cal fluctuations, as stated by the standard error of the mean, dominate the 
uncertainty of the hyperfine splitting {i>hf) measurement as seen in Table [H 

Table [2] presents the results of the hyperfine splitting of the 6Si/ 2 level as 
well as the corresponding values of the magnetic dipole constants for both 
isotopes. Our results are consistent with the past measurements of Gupta et 
al. [IT] and represent a precision improvement by a factor of 63 for 87 Rb and 
by a factor of 30 for 85 Rb. The table contains the theoretical predictions for 
85 Rb of Safronova et al. [16]. The theoretical calculation includes corrections 
for the finite size of the nuclear magnetic moment distribution, which is mod- 
elled as a uniformly magnetized ball. The theory value for 87 Rb comes from 
multiplying the theoretical value of the splitting in 85 Rb by the ratio of nuclear 
gyromagnetic ratios and does not include a hyperfine anomaly difference. 

Optical pumping effects are the most delicate of all the systematic effects. The 
polarization of the lasers determine the relative sizes of the peaks (m sublevels) 
that form the resonances of the 6S1/2 hyperfine levels. For our experimental 
conditions (around 1 G with circularly polarized light) we observe no difference 
between positive and negative magnetic field directions. Table [T] shows the 
bound for this effect. 



Systematic effects 


u HF ( 85 Rh) ( MHz ) 


u HF ( 87 Rb) (MHz) 


Optical pumping effects 


< 0.016 


< 0.029 


Power of 795 nm laser 


< 0.020 


< 0.005 


Power of 1.3 ^m laser 


< 0.011 


< 0.011 


Atomic density 


< 0.020 


< 0.010 


Non linear fit 


< 0.028 


< 0.023 


B-field fluctuations 


< 0.015 


< 0.025 


Total Systematic 


< 0.047 


< 0.047 


Statistical error 


0.100 


0.160 


TOTAL 


0.110 


0.167 



Table 1 

Error budget for the hyperfine splitting measurement 



6 



Rb [MHz] 87 Rb [MHz] 



vjjf this experiment 


717.540(110) 


1615.320(167) 


A previous experiment [T7] 


239.3(12) 


809.1(50) 


A this experiment 


239.18(4) 


807.66(8) 


A theory |16j 


238.2 


807.3 



Table 2 

Hyperfine splittings and magnetic dipole constants for the 6S1/2 level. 

We change the power of the 795 nm and 1.3 /zm lasers to look for any kind 
of dependence that would indicate line splitting effects such as the Autler- 
Townes effect or power broadening. We observe no systematic variation in 
the hyperfine separation as a function of laser powers. No systematic effect 
is found in the hyperfine splitting for changes in the temperature of the cell 
from 23 to 40 °C. The bounds for these effects are on Table [U 

The observed linewidths range from 30 to 40 MHz and fit (reduced y 2 ) Lorentzian 
profiles better (« 1) than Gaussian (« 10), with the Voigt profile between the 
two. We do not observe changes in the splitting that depend on the frequency 
range fitted around the resonances. The values obtained for the difference be- 
tween the centers of the resonances using Gaussian fits give consistent results 
with the Lorenzian fits, not changing the measured splitting by more than the 
reported error. 

The three layers of magnetic shielding limit the external magnetic fluctuations 
to less than 1 mG, while the stability of the magnetic field from the solenoid 
is comparable. We operate within the low magnetic field regime (1 G) and 
extrapolate to B=0 for both polarization sequences with linear fits with an 
average correlation coefficient of 0.97. 

Figure [3] shows the normalized ratio of hyperfine constants from Eq. [2] to show 
the hyperfine anomaly difference extracted from our measurement in the 6S1/2 
level of 85 Rb and 87 Rb as well as the corresponding values for other levels cur- 
rently in the literature. We use the values from Duong et al. [T5] to obtain 
the ratio of g values g 85 / 'g 8 7 = 0.295061 ± 0.00003, which is consistent with 
the values quoted by Stone [IB] , to calculate the hyperfine anomaly differ- 
ence from the measured A coefficients. The hyperfine splittings of the 55*1/2 
ground state are very well known as they are used in atomic clocks [19j . The 
5Pi/2 ratio comes from the work of Barwood et al. [20] which is consistent 
with more recent frequency comb measurements [2"Tf2"2"] . The 6S1/2 ratio cor- 
responds to the present measurement which gives 87^85 (6S1/2) = —0.0036(2). 
The difference in the anomalies is indeed a factor of thirty larger than the 
expected BCRS contribution and it comes from the BW effect. The 6P1/2 ra- 
tio comes from the critical evaluation of Arimondo et al. |19j . while we find 
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1.001 - 




0.996 " 



Fig. 3. Ratio of hyperfine constants normalized by the nuclear g factors showing 
hyperfine anomaly differences in 85 Rb and 87 Rb based on five different electronic 
states. The value for the §S\/2 (circle) comes from the present measurement. See 
the text for the other references (squares). 

from the recent frequency comb measurement of the 7S±/2 state by Chui et 
al. [23] 87^85 (7Si/2) = —0.0033(2), although they do not calculate nor mention 
the hyperfine anomaly in their paper. We have not been able to find in the 
literature values for higher levels with adequate precision to include them in 
the figure. Fig. [3] shows that the hyperfine anomaly difference measured with 
the 5*1/2 states is independent of the principal quantum number, just as Bohr 
and Weisskopf predicted pQ. The measurements are consistent with each other 
as are the P\j2 states that show a much smaller anomaly. 

We have measured the hyperfine splittings of the 6Si/ 2 level of 85 Rb and 87 Rb 
to a precision of 153 ppm and 103 ppm, respectively. Our measurement allows 
us to extract the hyperfine anomaly difference of the two isotopes with a 
precision of 5%. It is in excellent agreement with that of the ground state, 
confirming that it is the angular momentum and not the principal quantum 
number that enters in the Bohr Weisskopf effect [lj . We show also that nuclear 
magnetization information can be extracted through careful probing of excited 
state hyperfine splittings. 

Work supported by NSF. E. G. acknowledges support from CONACYT. 
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